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'~~7- Abstract 

<. 

*^ , Hypersingular integrals of the type 

t^ : UTn,m,r) = j'^ Tll^^^^r:ids , |r| < 1 



(s — rY 



and 



00 ! Ia{Un,m,r) = 4 — ds , \r\ < 1 



1 TT t,.\(A „2\m— 



>\ T nr ^ T Un{s){l-s 



1 



{s 



\^ • are investigated for general integers a (positive) and m (non-negative), where Tn{s) and 

^P . Un{s) are the Tchebyshev polynomials of the 1st and 2nd kinds, respectively. Exact formu- 

las are derived for the cases a = 1, 2, 3, 4 and ttt, = 0, 1, 2, 3; most of them corresponding to 
new solutions derived in this paper. Moreover, a systematic approach for evaluating these 
integrals when a > 4 and m > 3 is provided. The integrals are also evaluated as \r\ > 1 
in order to calculate stress intensity factors (SIFs). Examples involving crack problems 
are given and discussed with emphasis on the linkage between mathematics and mechan- 
ics of fracture. The examples include classical linear elastic fracture mechanics (LEFM), 
K> , functionally graded materials (FGM), and gradient elasticity theory. An appendix, with 

5h ' closed form solutions for a broad class of integrals, supplements the paper. 

1 Introduction 

Finite and boundary element methods are two of the most frequently used numerical approaches 
for solving crack problems in fracture mechanics. An alternative approach is the integral equa- 
tion method, which is more efficient {e.g. it reduces a partial differential equation (PDE) in two 
dimensions to an one- dimensional integral equation) and, in general, is more accurate than the 
aforementioned methods. The accuracy of the integral equation method relies on the analytical 
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evaluation of singular kernels to cancel the singularity (regularization) . In general, the cancel- 
lation of singularity is not trivial, in particular, in the case of hypersingular integrals. This is 
the main concern of this paper. 

Integral equations arising in static crack problems in fracture mechanics are typically Fred- 
holm integral equations of the form 

I kernel(x,t) D{t) dt = p{x) , c < x < d , (1) 

J c 

where kernel (x, t) is, in general, a singular function of (x, t); D{t) is the unknown, called density 
function; p{x) is some known (input) function corresponding to the loading on the crack faces; 
and the interval (c, d) refers to the crack surfaces where 2a = d — c denotes the crack length. 
By the Fourier transform, we write 

/oo 
K(Oe'(*-")«t/e • (2) 

■oo 

The singular part of the kernel can be separated from the regular part, by decomposing the 
Fourier transform as 

K{0 = £oo(0 + F(0 - /^oo(0] , (3) 

singular nonsingular 

which can be accomplished through asymptotic analysis (discussed later in this paper). Such 
an analysis is difficult for complicated K{C,). This is another issue to be addressed in this paper. 
Once the decomposition (|^) is accomplished, the integral equation (|l|) can be rewritten as 

T r — ^ dt + k{x,t) D{t) dt + /(x) = p{x) , c<x<d , (4) 

J c \t ^) J c 

where ^ denotes an improper integral; Ca is a constant associated to the singular kernel 
l/(t— x)"; k{x, t) is the nonsingular (regular) kernel; /(x) is a function standing for the free term; 
and a is a positive integer which determines the degree of the singularity. If a = 1, the integral 
equation (^ is called a Cauchy singular Fredholm integral equation, and the singular term is 
evaluated as a Cauchy principal-value (CPV) integral. If a > 2, it is called a hypersingular 
Fredholm integral equation and the singular term is evaluated as a Hadamard finite-part (HFP) 
integral [?, ?, ?, ?]. The notation-^ and -^ refer to CPV and HFP, respectively. Most of 
the works in the literature involve either a = 1 [?,?,?,?,?, ?] or a = 2 [?, ?, ?]. Thus 
another focus of this paper is to deal with hypersingular integral equations with a > 3 which 
arise naturally in gradient elasticity theories (see Example 3 in Section |7|). 

The singular and hypersingular integrals which involve Tchebyshev polynomials (T„, first 
kind; Un, second kind) and weight function (1 — s^)™'~2, -m > 0, with singularity a > 1 are 
defined by 

Ic,{Tn,Tn,r) = =h ds , |r| < 1 , (5) 



and 



Ia{Un,m,r) 



' Unis)il 



-.2\m- 



-ds 



\r\ < 1 



(6) 



The scope of this paper is as follows. First, Cauchy singular integrals, i.e. a = 1, are evaluated 
and exact formulas are derived for general m. The new results here are closed form analytical 
solutions for /i(T„, m, r), m > 1 and /i(t/„, m,r), m > 2. Once /i(T„, m, r) and Ii{Un, m, r) are 
known, hypersingular integrals Ia(Tn, m, r) and Ia{Un, m,r), a > 2, can be found by successive 
differentiation (with respect to r) in the sense of finite-part integrals; formulas for /2(T„,m,r), 
m > 1 and l2{Un, m,r), m >2 are derived in this manner. In the cases where a > 3, evaluation 
of hypersingular integrals becomes tedious and the formulas are lengthy. Thus Ia{Tn, rn, r) and 
Ia{Un, m, r) are provided only for a = 3, 4 and general m. 

2 Related Work 

Singular integral equations have played an active role in the field of solid mechanics, particularly 
in the solution of fracture mechanics problems. For instance, the Journal Integral Equations 
and Operator Theory has dedicated a special issue to "Integral Equations Methods in Engi- 
neering and Physics" (Volume 5, No. 4, 1982). Also, in June 1984, IMACS (International 
Association for Mathematics and Computers in Simulation) has held a Symposium [?] devoted 
to "Numerical Solution of Singular Integral Equations" . 



According to the notation introduced in equation (^), singular integral equations can be 
classified by the order of singularity a. The case a = 1 has been widely used and well devel- 
oped [?,?]. A rich field of application of singular integral equations is fracture mechanics of 
bimaterial and nonhomogeneous materials. For instance, the investigation of crack behavior 
in nonhomogeneous materials has found many applications to functionally graded materials 
(FGMs) [?, ?, ?]. Another use of singular integral equations involves FGMs for high tempera- 
ture applications, so that thermal stress intensity factors can be numerically calculated [?, ?]. 
Application of hj^ersingular Fredholm integral equations for a > 2 can be found in refer- 



ences 



r? ? ?i 



Quadrature formulas which involve hypersingular integrals have been drawing a considerable 
amount of concentration [?, ?, ?, ?, ?, ?, ?, ?] after Kutt first introduced the Hadamard finite- 
part (HFP) idea in his work [?]. Based on some previous work, Kaya [?] has presented a very 
nice interpretation about HFP integrals. One of the key steps in the derivations involves the fact 
that higher order singular integrals can be obtained from lower order ones by exchangeability 
of differentiation and integration [?, ?], e.g. 



D(s] 



-1 



\a+l 



-ds 



1 
a 



' I. 

I dr 



Dis) 



ds 



1 d 
a dr 



Dis) 



-ds 



-1 



where D{s) is the normalized density function. For instance, in order to find 

^' D{s) 



-1 (s-r)' 



;ds , |r| < 1 



\r\ < 1, 



(7) 



it suffices to know how to evaluate 

^' D{s) 



ds , |r| < 1 . 
I s — r 

This concept is apphed later in this paper. 

Another main motivation for numerical evaluations of hypersingular integrals is due to the 
boundary element method, and the reader is directed to the review paper by Tanaka et al [?]. 
Most recent work has been focused on singularity with a = 2 in two-dimensional problems. 

3 Theoretical Aspects 

First, relevant concepts involving integration and approximation are given. These concepts po- 
sition the contribution of the work with respect to the available literature. Next, a discussion 
on the influence of the density function on the corresponding singular integral equation formu- 
lation is presented. Afterwards, basic properties of the Tchebyshev polynomials are provided. 
These properties are heavily used in the analytical derivations that follow. 

3.1 Integration and Approximation 

As far as the integration and numerical procedures are concerned, the integral equation (^) may 
be normalized through the following change of variables 



d-c \ 1 ) d-c\ 1 _ 

which leads to the normalized version of the integral equation (Q) written asQ 



/ 7^^^^^+ \ }Cir,s)Dis)ds + F{r) = P( 
J-i [^ ^ ''") J -I 



-1 < r < 1 . (9) 



The density function D{s) is further assumed to have the representation 

D{s) = R{s)W{s) . (10) 

The weight function W{s) determines the singular behavior of the solution D{s) and has the 
form 

W{s) = (l-s)'"i(l + s)'"2 . (11) 

In general, mi ^ 1712, and the corresponding integrals, which involve Jacobi polynomials 

Pi'"^''"^^(s), areofthetype 

I s — r 



^The notations in this paper have been chosen as foUowing: x and t refer to the physical quantities and have 
dimension of "Length" ; r and s are normahzed (dimensionless) variables, corresponding to x and i, respectively. 



and can be expressed in terms of gamma and hypergeometric functions [?, ?, ?, ?]. In this 
paper, only the case nii = m2 is considered and mi, m2 are set to be 

mi = 1712 = Tn . (13) 

Thus W{s) can be expressed as 

W{s) = (l-s2)"'"5 m = 0,l,2,--- . (14) 



According to function-theoretic method [?, ?, ?, ?, ?] , the value of m is determined by the 
order of singularity a. For example, if a = 1, then m = and the fundamental solution D{s) 
to the Cauchy singular integral equation (^ takes the form 

D{s) = R{s){l-sr'^ = -|kL. (15) 

In this case, which consists of the majority of the work involving applications of integral equa- 
tions to fracture mechanics [?, ?, ?, ?, ?, ?] , R{s) is chosen to be 

OD 

R{s) = 5^a„T„(s) ; (16) 

n 

and because of that, the CPV integral /i(T„, 0, r) can be evaluated exactly [?,?]: 

/i(T„,0,r) = / ^f'l ds = 7r^„_i(r) , n>l. (17) 

J_i (s -r)Vl - s^ 

Another reason for choosing the approximation (p!5| ) is that with respect to the weight function 
W{s) = l/Vl — s^, the class of the Tchebyshev polynomials of first kind T„(s) is an orthogonal 
family [?, ?]: 

^ T ( ^T ( \ { ^ m = n = 

im {S)in{S ) ^^^ I ^^2 m = n; m,n = 1, 2, 3, ■ ■ ■ (18) 

1 vl-s^ [o m^n; m,n = 0, 1, 2, ■ ■ ■ 

With this orthogonal property a Galerkin-type method [?] may be applied to find the coefficients 
Qn in equation (0). 

If a = 2, then m = 1, and the solution D{s) to the hypersingular integral equation (|]) is 
characterized by 



D{s) = R{s){l-s^)2 = R{s)VT^^ . (19) 

Correspondingly, R{s) is chosen to be 

oo 
R{S) = Y,bnUn{s) , (20) 



because of the same reasons for the case a = 1, namely, analytical evaluation and orthogonal 
property. With respect to the first reason, the HFP integral /2(f/„,l,r) can be evaluated 
analytically [?] : 

/■^ U (s)\/l - s2 
/2(t/„,l,r) = 4 "^^'^ ,, ds = -in + l)nU^ir), n>0. (21) 

J-i [s-r)^ 

According to the second reason, by orthogonality. 



7r/2 m = n; m,n = 0,1,2, 



_^ '"v; ™v;v |Q m^n; m,n = 0, 1, 2, ■ • • , ^ ' 

and one may apply Galerkin-type methods [?] to find the coefficients 6„ in equation (pO|). 

When 771 = 3, then Vr(s) = (1 — s^)^/^, and neither T„(s) nor f/„(s) is an orthogonal family. 
However, if collocation method is applied, one does not need the orthogonal property, as long as 
the expansion function R{s) is chosen such that 



^ i?(s)(l-s^)5 



2^. 

ds 
_i (s-r)" 



can he evaluated analytically. For example, if R{s) is expandeded as a Tchebyshev polynomial 
of the 1st kind T„(s) or the 2nd kind t/„(s), i.e. 

oo oo 

R{s) = Y,^nTn{s) or R{s) = Y,bnUn{s) , (23) 

n n 

then the evaluation of 

T (T ^ f T„(S) (1 - sT~'^ . ,,,, , / [/„(.) (l-.T"^ . 

Ia{Tn,m,r) = f — ds or I^[Un,m,r) = 4^ — ds 

J -I (s-r)" J_i (s-r)" 

for general m = 0, 1, 2, ■ ■ ■ and a = 1, 2, 3, ■ ■ ■ is a necessary step for the numerical approach 
to the integral equation (|^). This is the one of main tasks in this paper and is addressed in 
Sections 4 and 5 . 

3.2 Selection of the Density Function 

Usually the unknown function D{t) in equation (|l]) can be chosen as the displacement pro- 
file {e.g. u{t) - a displacement function), the (first) derivative of the displacement function 
{du(t)/dt, denoted by (pit) - the slope function), or a higher derivative of u(t). The choice of 
the unknown function D{t) will affect the degree of singularity in the formulation. For example, 
consider the standard mode III crack problem in a free space [?] and a linear elastic fracture 
mechanics (LEFM) setting. If D{t) is chosen to be the slope function 0(t), then the governing 
integral equation is the Cauchy singular integral equation 

D{t) = <p{t) , 4 —^dt = p{x) , c<x<d . (24) 



However, if D{t) is chosen to be the displacement function w{t), then the hypersingular integral 
equation with a = 2 is obtained, 

D{t) = w{t) , 4 ^^ ' dt = p{x) , c<x <d . (25) 

The differences between the above two formulations are discussed next. 

In general, the (numerical) solution of a Cauchy singular integral equation, e.g. equation 
(^), is easier than the corresponding hypersingular integral equation, e.g. equation (p5|). A 
quick observation is that if equation ( P^ is used, then the actual crack surface displacements 
are obtained directly; while if equation (p5|) is used, an extra step of integration is needed to 
recover w{t). However, integration is not an unpleasant thing to do. 

So, what can be gained from a more singular equation such as (|25|) ? A formulation with 
more singular integral may lead to a simpler kernel function and, thus, simplify the kernel eval- 
uation and decomposition described in equation (|). Issues regarding the differences between 
formulations of the type given in equations ( pi] ) and ( P^D will be discussed in detail in the 
example section of this paper. 

3.3 Properties of Tchebyshev Polynomials 

The evaluation of Cauchy singular and hypersingular integrals which involve the Tchebyshev 
polynomials T„(s) and Un{s) highly depends on the special properties of these polynomials. 
They are listed here for the sake of completeness and because they will be of much use later in 
the development of this work. Most of them (but not all) can be found in Hochstrasser [?] and 
Kaya and Erdogan [?]. 

• Definition of Tchebyshev polynomials of the first kind: 

T„(s) = cos[?2cos~"'^(s)] , n = 0, 1, 2, • • • (26) 

• Definition of Tchebyshev polynomials of the second kind: 

^ sm|(r, + l)cos--W| _ „^„_j_2,... (27) 

sm COS ^ s 



Iterative (recursive) properties: 

~ 2 



sT„(s) = i[T,,+i(s) + T„_i(s)] , n > 1 (2^ 



sUn{s) = ^[Un+i{s) + Un-iis)] , n > 1 (29) 



T„(s) = i[f/„(s)-t/„_2(s)] , n>2 (30) 



f/„(s)(l - s^) = sT„+i(s) - T„+2(s) , n>0 



(31) 



By means of equation (|28D, one may rewrite equation (31) above as 

Un{s) = ^^^^^[Tn{s) - T^+2{s)] , n > 



(32) 



Thus an additional equality, which is useful in handling cubic hypersingular integrals can 
be derived^: 

(32) 1 



u,,is)ii - s')^ ''=' - [r„,(.) - r„+2(5)] VT^^ 



-f/n+2(s) - f/„(s) + -f/„_2(s) 



Vl-s2 , n> 2 



[f/„+2(s) - 2Un{s) + f/„_2(s)] Vl-s2 , n > 2 



(33) 



Derivatives: 



dTJs) 



ds 



nUn-i{s) , n > 1 



cff4(s) 
ds 



1-S2 



— ^f/„_i(s) - -f/„+i(s) 



ra > 1 



(34) 
(35) 



4 Cauchy Singular Integral Formulas {a = 1) 

This section mainly evaluates Ii{Tn,m,r) and Ii{Un-,Tn,r), which are defined in equations (^ 
and (P). The new result here is that the singular integral formulas are found for general 
m. In order to obtain this new result, two well known Cauchy singular integral formulas are 
introduced [?, ?]: one is already stated in equation (p!?!), and the other one is 



/i(t/„, 1, r) = f ds = -7rT„+i(r) , n > 

J -I s-r 



(36) 



which can be obtained as follows 



/i(f/„,l,r) = f ds 

^1 s — r 



^ 1 [^ Tn.(s) -T„.+2(s) 

Vl — s^(s — r) 
^ ^|f/n-,(r)-(7„+,(r)l 

® ->rT,.„(r). 



ds 



^The equation number is stacked above the equal sign to show how the equations are being derived and 
connected. 



The integral formulas for ?7i = 0, 1, 2, 3, ■ ■ ■ are derived below. The general formulas have the 
restriction of minimum n. The lower n terms can not be derived by general formulas, and are 
given in Appendix A. 

4.1 Ii{Tn,m,r), m = 0,1,2,3 

• /i(T„,0,r): 

This is equation (|1^. 



/i(T„,l,r) 



' T^{S)VT^^^^ (30) 1 / [^n(^)-^n-2(^)]v^r^^ ^„ 



_i s — r 2 J_i s — r 

® |[r„_i(r)-r„+i(r)], n>2. (37) 



/i(T„,2,r) 



' T„.(.)(l-.^)f ^^^ (30) f llUnjs) - Un-2{S)]{1 - S')-^ ^^ 



1 s — r j_i s — r 



® I {[T„_i(r) - 2T„+i(r) + T^+,{t)] - [T„_3(r) - 2T„_i(r) + T„+i(r)]} 



4[T„_3(r) - 3T„_i(r) + 3T„+i(r) - T„+3(r)] , n > 4 (38) 



/i(T„,3,r): 



1 T„(s)(l-g^)l ^^^ (30) / |[[/.(s)-[/„_2(3)](l-g^)^ ^„ 



_i s — r j_i s — r 

^ ^ {[T„_5(r) - 4r„„3(r) + 6T„_i(r) - 4T„+i(r) + T„+3(r)] - 

[r„_3(r) - 4r„_i(r) + 6T„+i(r) - 4T„,+3(r) + T„+5(r)]} , n > 6 
= ^[Tn-rAr] - 5T„_3(r) + 10T„_i(r) 

-10T„+i(r) + 5T„+3(r) - T„+5(r)] , n > 6 (39) 



4.2 Ii{Un,m,r), m = 1,2,3 

• /i(t/n,l,r): 
This is equation 

• /i(f/„,2,r): 



1 f/„(s)(l-s2)l^^(32) 1 / [T„(s)-T„+2(g)]v^r^^ ^^^ 



1 s — r 2 y_;^ s — r 

® ^ {[T„_i(r) - T^+,{t)] - [Tr^^^ir) - T^^r)]} , n>2 
= ^[T„_i(r) - 2r„+i(r) + r„+3(r)] , n > 2 (40) 



/i(f/n,3,r): 



' Unis)il-S')l (32)1 f [T„(s)-T„+2(S)](1-S')i 



-ds 



ds 



I s — r 

g) ^ 
16 

[T„_3(r) - 3T„_i(r) + 3T„+i(r) - 3T„+3(r)]} , n > 4 



{[T„_i(r) - 3T„,+i(r) + 3T„+3(r) - T„+5(r)] - 



vr 



---[T„_3(r) - 4T„_i(r) + 6T„+i(r) - 4T„+3(r) + T„+5(r)] , n > 4 (41) 
lb 



4.3 /i(r„, m,r) and /i([/„,m, r) 

At this point one may easily see the procedural steps above, which take advantage of recursive 
properties ( pOD and ( |32D between the Tchebyshev polynomials T„(s) and Un{s). For instance, 
evaluation of Ji(T„,4, r) = j^-^ T„(s)(l — s'^y/'^/{s — r)ds can be reduced to evaluation of 
/i(f/„,4, r) = ^_^ Un{s){\ — s^)^/^/(s — r)(is, which, in turn, can be reduced to evaluation of 

/i(T„,3,r) = ^_^ T„(s)(l — s^)^/^/(s — r)(is. After a suitable number of steps, this reduction 
leads to either ([T7|) or (|36|) . This procedure is summarized in Figure |l|. 
Thus, by induction, one obtains the following formulas for ?7i > 1. 

• Ii{Tn,m,r), where m>l, and n > 2m 



J-1 s-r 


= 7r(-l)™+^ 


.,x 2m-l2m-l 


/2m-lN 


Tn+l-2m+2j{r) ■ 










(4 



2) 



Ii{Un, m, r), where m > 2, and n > 2m — 2 



f/n(s)(l-s2) 



2\m- 



-(is = vrf— 11 



2m-2 2m--2 



EM 

j=0 



2m-2 
*^ ' ^ Tn+3-2m+2i('^) • 



(43) 



The usual notation 



[m] 



]\ [m- ]^ 



denotes the binomial coefficients. 



5 Hypersingular Integral Formulas (a > 2) 

Once a Cauchy singular integral formula has been reached, all other hypersingular integral 
formulas may be obtained successively by taking differentiation with respect to r, and making 
use of the finite-part integral formula (|^). 



10 



rl Tn(g)(l-g2)"^+2 
J—1 s—r 

J| (H) 

1 T^(5)(i-52r-2 



(i5 



ds 



f: 



1 



5— r 



ds 



fl. wiizfriirf. 



1 



s—r 

J| (H) 



Equation ( [tt] ) or (gQ) 



Figure 1: Evaluation of Ii(Tn,m,r) and Ii{Un,m,r) for general m. The procedure reduces 
integrals with higher m to those with lower m. 



5.1 l2{Tn,m,r) 
By means of 



' Tn{s){l-S 



2\m- 



-ds 



d f' T^{s){i - s'r--2 



-1 



[s — r] 



dr 



ds 



s — r 



one readily obtains: 

. /2(T„,0,r) [?]: 



Tnis) 



1 VT 

(35) vr 



-(is = vr 



s^ s — r 



dUn-ijr) 
dr 



l-r2 






n>2 



(44) 



11 



/2(T„,l,r): 



as = - — [Tn-i{r) - T„+i(r)J 



1 

(34) vr 



s — r 



2dr' 



[(n - l)f/„_2(r) - (ri + l)f/„(r)] , n>2 



(45) 



HT„,2,r) 



I (s — r)2 

(34) TT 



ds = -^-rlTn-sir) - 3T„_i(r) + 3r„+i(r) - T„+3(r)] 
8 ar 



[ (n - 3)t/„_4(r) - 3(71 - l)f/„_2(r) + 3(n + l)?7„(r) 
(n + 3)f/„+2(r) ] , n > 4 



(46) 



/2(T„,3,r): 



-ds 



-1 



s — r 



= ^^[^-"5(r) - 5T„_3(r) + 10T„_i(r) - 10T„+i(r) + 5T„+3(r) - T„+5(r)] 
=^ ^ [(n - 5)t/„_6(r) - 5(n - 3)f/„_4(r) + 10(n - l)f/„„2(r) - 10(n + l)f/„(r) 



+5(n + 3)f/„+2(r) - (n + 5)f/„+4(r) ] , n > 6 



(47) 



l2{Tn, m, r), where m > 1, and n > 2m + 1: 



f-\ 



1 T^(g)(l-g2)"^-V2 ^ 

^ «-5 



s— r 



7r(-l)-+i(i)^"-^Vto (-1)^' 



2m- 1 



(n + 1 - 2m + 2j)Un-2m+2j{r) 



(4^ 



5.2 l2{Un,m,r) 
The following equality 



1 ?7„(.)(i-3^)-i ^^^ = 1 i' ^»(s)(i-s^r"^ ^„ 



s — r 



c^'" j-i 



s — r 



leads to: 



' ^"W^^^d. = -.^IM^a-.(„ + i)c/.M, „>0 



-1 



s — r 



dr 



which is the same as (W^ 



(49) 
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/2(f/n,2,r): 



' Un{s)il-s^)l 



n 



d 



[s — r] 



ds = -j-rl'^n-iir) - 2T„+i(r) + T„+3(r) 
4 ar 



(34) TT 



(n - l)f/„-2(r) - 2(n + l)f/„(r) + (n + 3)f/„+2(r)] , n > 2 (50) 



/2(f/n,3,r): 

^1 Un{s)il-S')l 



ds 



1 (s - ry 

= -T7^:r[^"-3(^) - 4T„-i(r) + 6T„+i(r) - 4T„+3(r) + T„+5(r)] 
lo ar 

('34'l TT 

= -- [{n - 3)f/„^4(r) - 4(n - l)f/„_2(r) + 6(n + l)f/„(r) - 4(n + 3)f/„+2(r) 
lo 



+ (n + 5)f/„+4(r) ] , n>4 



(51) 



l2{Uni fn, f^), where m >2, and n > 2m — 1: 



f-\ 



1 u^{s){i-s^r-i , 

7 T^ as 



TT 



2m- 2 



-ir (I)"""' E?ro"'(-l)^' ( , ) (n + 3 - 2m + 2j)f/n+2-2m+2,(r) 



(52) 



5.3 h{Tn,m,r] 
By means of 



' Tn{s){l-S 



2\m- 



1 d f T„(s)(l-s 



-1 



one obtains: 



(s — r)^ 



-(is = --— , 
2 dr y_i 



s — r 



2\m-i 

^ ds ^ 



h{Tn,0,ry. 



' T„(s)v/r^7 



ds 



TT 



IS — r] 



(^ + l)(n + 2)f/„_3(r) 



.(1_^2)2 

-2{n^ - 3)Un-i{r) + (n - l)'t/„+i(r)] , n > 3 



(53) 



/3(T„,l,r): 



^ Tn{s)VT^^ 



ds 



TT 



[s — r] 



:i-r2 



(n^ -n)f/„„3(r) 



(2n2 + 2)f/„_i(r) + {n^ + n)f/„+i(r)] , n > 3 



(54) 
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ds = ^^r^A -{n + S){n + 2)U^^,{r) 



-1 (s-r)3 "" 32( 

+ [{n + 3)(n + 4) + 3n{n + l)]f/n+i(r) - [3{n + l)(n + 2) + 3(n - l)(n - 2)]f/„ 

+ [3n(n-l) + (n-3)(n-4)]f/„_3(r)-(n-3)(n-2)f/„_5(r)} , n>5 



(55) 



J3(T„,3,r): 



ds = ^^^J_ [in' + 9n + 20)f/„+5(r) 



_i (s-rf '' 128(1- 

- 6(n' + 6n + 10)[/„+3(r) + 15 (n' + 3n + 4)f/„+i(r) 

- 20(^2 + 2)Un-i{r) + 15(n' - 3n + 4)f/„_3(r) - 6(n' - 6n + 10)f/„_5(r) 
+ (n^ - 9n + 2)f/„_7(r)] , n>7 



IsiTn, m, r), where m > 1, and n > 2m + 2: 



#-!. 



1 Tn{s){l-S^r-^ , 
? 71 "S = 



'2m-l' 

-1)M ) (n+l-2m + 2j)x 



[(n + 2 - 2m + 2j)f/„_i_2m+2j(r) - (n - 2m + 2j)f/„+i_2m+2i(r) 



5.4 h{Un,m,r] 

By means of 



1 [/„(.)(i-5^r-^ ^_irf / u^{s){i-s'r~'2 ^^ 

I (s — r)^ 2drJ_i {s — r^ 



one gets: 

• /3(f/n,l,r): 



1 Un{s)VT^^ 



ds 



-1 [s-r) 

TT 



4(1 -r2) 



hiUn,2,ry. 



' Un{s){l-S')l 

as 



{n' + 5n + 6)Un+3{r) 



1 (s-r)^ 16(1 -r2) 

+ (3n^ + 9n + 12)f/„+i(r) - (3n^ + 3n + 6)f/„_i(r) 
+ (n^ - n)Un-3{r)] , n>3 



(56) 



(57) 



{-{2n' + 3n + 2)f/„_i(r) + (n" + n)Un+i{r)} , n > 1 (58) 



(59) 
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/3(f/n,3,r): 



ds 



n 



[s — r) 



n^ + 9n + 20)f/„+5(r) 



64(1 -r2) 

- (5n2 + 31n + 54)f/„+3(r) + (lOn^ + 34n + 48)f/„+i(r) - (lOn' + 6n + 20)f/„_i(r) 
+ (5n2-lln+12)f/„_3(r)-(n2-5n + 6)f/„_5(r)] , n>5 (60) 

l3{Un,Tn,r), where m>2, and n > 2m: 



^-!. 



1 Un{s){i-s^rA, 

^ ds 



[s—r] 



,2m 



-ir(i)^™T^E?:o"'(-iv' 



2m- 2 



(n + 3-2m + 2j)x 



[(n + 4 - 2m + 2j)Un+i-2m+2j{r) -{n + 2-2m + 2j)Un+3-2m+2j{r)] 



5.5 U{Tn,m,r] 

By means of 



1 T„(s)(l-s2r-2^^ Id f T„(s)(l-s2)™-^ 



'-1 (s-r) 

one reaches the following results: 
• /4(T„,0,r): 



TJs) 



ds 
4 3rfr 



-1 



(s — r)^ 



(is 



1 [s — r l"*Vl — s 



-ds 



n 



■ [{n^ + 6n^ + lln + 6)f/„_4(r) 



4,/i _ .2--- 48(1 -r2)3 
-{3n^ + 6n2 - 25n - 44)f/„_2(r) + (3n^ - 5n^ - 19n + 37)f/„(r) 
-(n^ - Sn^ + 7n - 3)f/„+2(r)] , n > 4 



/4(T„,l,r): 



(is 



TT 



1 (s_r)4 "" 48(1 -r2)2 

(n^ - n)f/„_4(r) - {3n^ + 9n + 12)f/„„2(r) + (3^^ + 9n - 12)f/„(r) 
-(n2-n)f/„+2(r)] , n>4 



h{Tn,2,ry. 

^' T„(s)(i-s2; 



-(is 



TT 



(n^ + 6n2 + lln + 6)f/„+4(r) 



1 (s-r)4 192(1 -r2)2 

- (Sn^ + 18^2 + 43n + 30)f/„+2(r) + (lOn^ + 12^^ + 134n - 36)f/„(r) 

- (lOn^ - 12^2 + 134n + 36)f/„_2(r) + (Sn^ - 1872^ + 43n - 30)f/„_4(r) 

- (n^-6n2 + lln-6)f/„_6(r)], n>6 



(61) 



(62) 



(63) 



(64) 
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ds 



[s — r] 



TT r 1 47 

384(1 _ ,2)2 [ - ( 2^' + 6^' + y ^ + 30)f/n+6(r) 



7 197 

' :n^ + 30n^ + — -n + 120)f/„+4(r) 



21 327 

,2-^ . -^- . 2 — ----— :yn3 + 54n2 + — n + 180)f/„+2( 

+ (^^3 ^ 30^2 + ^n + 90)f/„(r) - (^n=^ - SOn^ + ^^ " 90)f/„_2(r) 

21 327 7 197 

+ (yn^ - 54n2 + — n - 180)f/„„4(r) - {-rv" - SOn^ + — n - 120)[/„_6( 
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+ (-n^-6n' + yn-30)f/„_8(r)J, n> 
Ii{Tn-, m, r), where m > 1, and n > 2m + 3: 



£i 



1 Tn{s){l-S 



2\Tn- 



\s—r] 



-ds 



.^Am+l (l\'i'"i+2 1 



(^)- 



E 



2m- 1, 



-IV 



2m- 1 



3 (l-r2)2 Z^j=0 

{[(n + 2 - 2m + 2j)(n + 3 - 2m + 2j)]f/„-2-2m+2,(r) 
[2(n - 2m + 2jf + 4(n - 2m + 2j) - 6]f/„_2m+2i(r) + 

[(n - 2m + 2j){n - 1 - 2m + 2j)]f/„+2-2m+2i(^)} 



(n + 1 - 2m + 2j)x 



5.6 l4{Un,m,r): 
By means of 



-1 



1 [/„(g)(l-g^)"^-^ ^_ 1 rf f' Un{s){l-S 

(s — r)"^ 3drJ_i (s — ry 



2\m- 



-ds 



one obtains 

• h{Un,l,r): 



ds 



IT 



,2\2 



^ u^{s){i-s'y2 

1 (s-r)4 "^ 24(1 -r^^ 

-(2n^ + 9n2 + lln + 6)Un-2{r) + (3n^ + 3n'^ - 2n - 6)f/„(r) 
-(n3-n)f/„+2(r)] , n>2 



h{Un,2,ry. 



' Un{s)il 



S ]2 



-ds 



n 



rr- + 6n^ + lln + 6)f/„+4(r) 



_i v^-<) 96(1 -r2)2 

{An^ + 18n^ + 44n + 30)f/„+2(r) + (6n^ + ISn^ + 54n + 42)f/„(r) 
(4n^ + Gn^ + 20n + 18)f/„_2(r) + (n^ - n)f/„_4(r)] , n > 4 



(65) 



(66) 



(67) 



(6^ 
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-1 



[s - ry 



-ds 



n 



47 



r 1 ^ ( 

192(1 -,-p[-<2"' + «"' + Y" + '2°'«»«'''' 



1 

2' 

+ (3^3 + 27^2 + 93n + 117)t/„+4(r) - (yn=^ + 45^^ + ^n + 165)f/„+2(r) 

15 105 

+ (lOn^ + 30n^ + UOn + 90)?7„(r) - (yn^ + — n)t/n_2(r) 

+ {3n^ - 9n^ + 2ln - 15)f/„_4(r) 

- (in3 + 3n2 + yn-3)f/„_6(r)], n>6 

liiUn, fn, r), where m >2, and n > 2m + 1: 



(69) 



i 



1 Unis)il-s') 

1 



2\m— 



1 



-(is 



2m~2 



(s-r)4 

(-ir(l)"""''l(T^E,to'^(-lV'^ ^. J(n + 3-2m + 2j)x 

{[(n + 4 - 2m + 2j){n + 5 - 2m + 2j)]f/n-2m+2i(r)- 
[2(n - 2m + 2j)^ + 10(n - 2m + 2j) + 10]f/„+2-2m+2i(r) + 

[(n + 2 - 2m + 2j)(ri - 1 - 2m + 2j)]Un+4-2m+2j (r)} 



(70) 



6 Evaluation of Stress Intensity Factors (SIFs) 

An important task is to evaluate the stress intensity factors (SIFs) at both crack tips, since 
the propagation of a crack starts around its tips. In mode III crack problems, SIFs can be 
calculated from 



Kin{d) = lim A/27r(x — d)ayz{x, 0) , [x > d) 



x^d+ 



and 



Kiii{c) = lim \/2'k{c — x)ayz{x, 0) , (x < c) 



(71) 



(72) 



Note that the limit is taken from outside of the crack surfaces and towards both tips. Usually 
the left hand side of integral equation (^) is the expression for cryz{x, 0) which is valid for x is 
inside the crack surfaces (c, d) as well as outside of (c, d). Thus to calculate SIFs, the key is to 
evaluate the following integrals which are obtained after proper normalization and the change 
of variables described in equation (|]) , 



and 



Sa{Tn,m,r) 



Sa{Un,m,r) 



1 T„(s)(l-s2)™-(l/2) 



s — r 



ds, r^(-l,l) 



' f/„,(s)(l-s2)--(l/2) 



ds 



r^(-l,l). 



(73) 



(74) 
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Note that the above integrals are not singular as x 7^ t for t G (c, d) and x ^ (c, d). 

The strategy to evaluate Sa(Tn,m,r) and SaiUn^rn^r) for general integers a (positive) and 
m (non-negative) is similar to the process for evaluating laiTn, m, r) and Ia{,Un, "o^, r). It consists 
of evaluating the integrals Si{Tn,m,r) and 5'i(f/„, m,r) by means of the reduction procedure 
described in Section 4.3, and taking differentiation (with respect to r) to obtain Sa{Tn,m,r) 
and Sa{Un, rn, r) for a > 2. The relevant derivations are provided below where the range of r 
is restricted to \r\ > 1. These formulas are used in calculating SIFs for the examples presented 
in Section 7 

6.1 5'i(r„,m,r) and 5'i(t/„,m,r) 

• S'i(T„,0,r) (This is a well known integral [?]): 



■ T„W ,^_^(r-vg5T|r|/r)"_ ^^^^ ^^^ 



1 (s — r)-\/l — s^ V^^ — 1 



^ ,j. 



S'i(Tn,l,r) 



'^iW^IHZ,, 



s — r 



(30) 1 f'Un{s)VT^^, 1 /■' f/„-2(s)v^-S 



(is / ^^ ^ " (is 



Si{To,2,r) 



2 y_;^ s — r 2 J_i s — r 

® TT^V^^^^T ("r - Mv^^^^l'j ^ ^ > 2 . (76) 



' Z£Mi^i!)irf, = ^(,2 _ 1) f^ _ My;:^3I ) . (77) 



^i(ri,2,r): 



^ ri(s)(l-s^)2 . vr 



s — r 2 V r 



rfs = -(r^ - 1) r - ^Vr2 - 1 . (78 



Si{Tn,2,ry. 



^T„(s)(l-s2)f^ iM3) 



vr r 



s — r 



rfs^^ M(^^2_^^|5 r-^Vr2-l ,n>2. (79) 
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Si(U,„l,r): 



•coVEZ,, 



1 



s — r 



m 1 r Us) .,3_ir w^^,^ 



Si{Un,2,r) 



r-i / I I \ "^"'" 

(B _7r(r-^V;^^^Tj , n>0. (80) 



1 s — r 

(32) 1 f'Tn{s)VT^^_^^_l f''In+2{sWT^^^ 



2 y_]^ s — r 2 y_i s — r 

n+l 



® 7r(r2-l) (r-^V^^^^T) , n>2. (81) 



The formulas for Si{Tn,m,r) and Si{Un-,Tn,r) with general ?7i can be deduced by the same 
procedure described in Figure |I|, and are listed below. 

J-i s-r 

(II \ " 

r-^VH^^l ,m>0andn>2m. (82) 

r / 



Q ,,, , r ^n(^)(l-S^)-^/^ 

J-i s-r 

(II \ " 

r _ ^Vr2 - 1 ] , m> 1 and n>2m-2 . (83) 

6.2 S2{Tn,m,r) and S2{Un,m,r) 

Differentiating (with respect to r) the formulas for 5*1 (T„, -m, r) and 5*1 (f/„, m, r), we obtain the 
formulas for 5'2(T„,m, r) and S2{Un,rn,r). 

- as 



1 



(s — r)^ 



n+l / I I \ n+l 

r\ 



= -7c{n + 1) — Vr2 - 1 ( r - — Vr^ - 1 J + 27rr I r - — Vr^ - 11 , n > (84) 

and 

/•^ T„(s)(l-s2)i , TTlrL 2 ,3/ \r\ r-. Y , . 

/ ^^^-^ '—ds = ^ r^-l 5 r-^Vr2-l ,n>2. 85 

y_]^ g _ ^ ^ \^ r / 
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6.3 S3{Tn,m,r) and Si{Un,m,r) 



The following formulas are obtained by differentiating twice (with respect to r) the correspond- 
ing formulas obtained in Subsection 6.1. 



^Un{s)^-S'Y 



■ds 



IT 

2 



s — r 

.2 



n' + 2n + 3)~ 3{n + 1] 



V^^^^I, 



r 



r Vr^ — 1 

r 



n+l 



n > . 



(86) 



' Ti{s){l 



s') 



-ds 



3tt 
T 



vr^ 



Vr 



(87) 



' To{s){l 



s2)f 



-1 



ds 



377 
T 



-v/^^^^T| ( 1 



v/f^^I 




n+1 



[n^ + 2n + 3)-3{n + l] 



— \ r Vr^ — 1 

r 



n-l 



(n2-2n + 3) - 3(n - 1) 



v/f^^T 



Vf^^I 



ra > 2 



(89) 



7 Examples 

Three examples are presented here, which emphasize various aspects of Fredholm singular 
integral equation formulations and their linkage to fracture mechanics. These examples are: 

1. Internal mode I crack in an infinite strip. 

2. Mode III crack problem in nonhomogeneous materials. 

3. Gradient elasticity theory applied to a mode III crack. 

The first and last examples consider homogeneous materials, and the second example con- 
siders nonhomogeneous materials, which has relevant applications to the field of functionally 
graded materials [?, ?]. The first two examples are from classical elasticity, and the last one 
is from gradient elasticity theory. The first example involves mode I cracks and the last two 
examples involve mode III cracks. All the examples are formulated by using hypersingular 
integral equations. For the first two examples the order of singularity a is 2, and for the last 
example a is 3. A detailed comparison between f/„ and T„ representations is given in the first 
example. A discussion on the influence of the density function on the order of singularity of 
the integral equation is presented in the second and third examples. The description of the 
examples is summarized in Table [^ 
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Table 1: Description of the examples. 




Description 


Example 1 


Example 2 


Example 3 


Homogeneous material 


V 




V 


Nonomogeneous material 




V 




Classical elasticity 


V 


V 




Gradient elasticity 






V 


Crack mode 


I 


III 


III 


Density function 


displacement (v) 


displacement (w) 


slope (0) 


Degree of singularity a 


2 


2 


3 


Weight function exponent 


m-(l/2) 


1/2 


1/2 


3/2 


Representation 


Un,Tn 


f/n, T„ 


Tn 



7.1 Internal Mode I Crack in an Infinite Strip [?] 

Consider a crack in an infinite strip of homogeneous material, as illustrated by Figure 
governing partial differential equations (PDEs) and boundary conditions are: 



The 



V2 / \ I 2 / d ulx,y) I a ulx,y) \ ^i 

uix,y) + — ( ^^ + ^iaf j = , 

V2 / \ I 2 / d^v(x,y) I d^v(x,y) \ n 

^v{x,y) + — i^-^^ + ^^j =0, 

crxxiO,y) = (7xyiO,y) = (^xxih,y) = a^yih,y) 



a 



xy 



(X, 0) = 



O-yy{X,0) = -P{X) 

vlx,0) = , 



— oo < X, y < oo , 

— oo < X, y < oo , 

— cxD < y < oo , 
< X < h , 
X G (c, d) , 

X ^ [c, d] , 



(90) 



where u and f are the x and y components of the displacement vector; ctjj is the stress tensor; 
K, is an elastic constant (k = 3 — 4z/ for plane strain, /t = (3 — z^)/(l + z^) for plane stress, and 
z/ is the Poisson's ratio.) This problem has been studied by Kaya and Erdogan [?] by means 
of a Un representation, and it has also been used as a benchmark problem by Kabir et. al. [?] . 
Here both Un and T„ are employed and compared. 



The governing integral equation can be written in the form given by equation @) as [7 



/\v{t) 



dt+ k{x,t)Av{t)dt = -IT 



1 + K 



p{x) , c < X < d 



(91) 



(t-x)2" ■ y, ■--'- ^-'"" •■ V 2/x 

where the primary variable is the crack opening displacement Av given by 

Av{x) = v{x, 0"*") — v{x, 0~) , c < X < d , 

and the kernel k{x,t) is given in Kaya and Erdogan [?], equations (51) - (54c), page 112. It is 
worth noting that as /i — >■ cxo (see Figure §), the integral equation for the half plane is recovered 
and the kernel k{x,t) is reduced to a much simpler formQ 

,, , -1 12a; 12x2 

k[x,t) 



{t + xy (t + a;)3 {t + xy ■ 



■^To be consistent with the notation adopted in this paper, we have used symbols different from reference [?]. 
For instance, in [?] upper case K{t,x) is used for k{t,x). 
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d 



X 



h 
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Figure 2: A mode I crack in an infinite strip. 

After normalization, tlie corresponding integral equation can be written in a fashion similar to 
equation (|^), i.e.Q 



c^ D(s) r^ 

+ - — ^-^ds+ }C(r,s)D(s)ds = P(r) , -1 < r < 1 , 
J_i [s-ry J^i 

where D{s) is the unknown displacement function, the regular kernel is 

, 12 [.+ (^)] 12 [.+ (^)]^ 



(92) 



/C(r, s) 



-1 



[(- + ^) + 2(^)]^ [(r + .) + 2(^)]^ [(r + s) + 2(ft.)] 
and the loading function is 



4 ' 



P(r) 



1 + K 

-vr I —i 1 p 



2fi 



d — c\ d + c 



The case c > represents an internal crack, which is the case of interest in this work. Based 
on the dominant behavior of the singular kernels of the integral equation (p2D, the solution 
takes the form 

D{s) = R{s)Vl - s2 . 



* Again, the notation is different from [?]. 
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1.6 



1.2 



+ 



0.8 



0.4 







- Y^- ^-^ ^j-^ - 

c = OA \d = 20.1 







10 



X 



15 



20 



25 



Figure 3: Displacement profiles for a mode I crack in an infinite strip obtained by means of f/„ 
and T„ representations (A^ + 1 = 15). Here c = 0.1, d = 20.1, 2a = 20, and {c+d)/{d—c) = 1.01. 
The crack is tilted to the left because of the "edge effect". 



Here the representation function R{s) is approximated in terms of Tchebyshev polynomials of 
1st and 2nd kinds, i.e. 



N 



N 



Ris) = ^a„f/„(s) and R{s) = ^^bnTnis) . 



n=0 



n=0 



The unknown coefficients a„ and 6„ are determined by selecting an appropriate set of collocation 
points 



Tj = COS 



(2n- l)7r 
2{N + 1) 



j = 1,2,- ■ ■ ,N + 1; for f/„ representation . 



r, = cos 



nil 

N + 2 



, j = 1,2,- ■ ■ ,N + 1; for T„ representation . 
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Once the solution is obtained, the SIFs can be calculated from^ 



Kj{c) = lim ^/2tt{c — x)ayy{x, 0) , (x < c) 

, (x > c) 



1 + k/ x^c+ ^27r(x - c) 



2/i 



d — c 

271 



i?(-l) 



(93) 



and 



Kj{d) = lim ^/2tc{x — d)ayy{x, 0) , (x > rf) 



2/i 



lim 



■yy\ 
D(x) 



1 + kJ x-.d- ^2Tx{d-x) 
2/i 



, {x < d) 



1 + K 



d — c 
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Ri+l) 



(94) 



which are obtained from equation (pT|) by observing that its left-hand-side gives the stress 
component ayy{x,0) outside the crack interval {c,d). 



Table 
terms 



2: Normalized stress intensity factors (SIFs) for an internal crack in a half-plane. A^ + 1 
are used in approximating the primary variable. 



d+c 
d-c 


A^+1 


Un Representation 


Tn Representation 


Kaya and Erdogan [?] 


Ki{c) 


Kj{d} 


Kj{c) 


Ki(d) 


Ki(c) 


Ki(d) 


poVT(d-c)/2 


P0y/7T{d-c)/2 


poy/7r{d^c)/2 


po-\/7r(d-c)/2 


po^Ti{d-c)/2 


P0y/TT{d-c)/2 


1.01 


15 


3.6437 


1.3292 


3.8037 


1.3313 


3.6387 


1.3298 


1.05 


10 


2.1541 


1.2535 


2.1920 


1.2543 


2.1547 


1.2536 


1.1 


10 


1.7583 


1.2108 


1.7655 


1.2111 


1.7587 


1.2108 


1.2 


6 


1.4637 


1.1625 


1.4728 


1.1632 


1.4637 


1.1626 


1.3 


6 


1.3316 


1.1331 


1.3346 


1.1335 


1.3316 


1.1331 


1.4 


6 


1.2544 


1.1123 


1.2556 


1.1125 


1.2544 


1.1123 


1.5 


4 


1.2036 


1.0966 


1.2066 


1.0969 


1.2035 


1.0967 


2.0 


4 


1.0913 


1.0539 


1.0916 


1.0540 


1.0913 


1.0539 


3.0 


4 


1.0345 


1.0246 


1.0346 


1.0246 


1.0345 


1.0246 


4.0 


4 


1.0182 


1.0141 


1.0182 


1.0141 


1.0182 


1.0141 


5.0 


4 


1.0112 


1.0092 


1.0112 


1.0092 


1.0112 


1.0092 


10.0 


4 


1.0026 


1.0024 


1.0026 


1.0024 


1.0026 


1.0024 


20.0 


4 


1.0006 


1.0006 


1.0006 


1.0006 


1.0006 


1.0006 



Table |^ presents the SIFs at both tips of an internal crack in a half-plane {h -^ cx)) under 
uniform load (p(x) = po) obtained with both [/„ and T„ representations. First, it is worth noting 
that the present SIF results for the Un representation compare well with those reported in Table 



^Kaya and Erdogan [?] do not consider the factor -0? in the definition of SIFs, equations (|92 
that this does not affect the normahzed SIFs {e.g. see Table g). 
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Figure 4: Displacement profiles for a mode I crack in an infinite strip obtained by means of [/„ 
and Tn representations (A^ + 1 = 8). Here c = 1, (i = 3, 2a = 2, and (c + d)/{d — c) = 2. 

1 (page 114) of the paper by Kaya and Erdogan [?] for the entire range of values describing 
the relative position of the crack, i.e. 1.01 < {d + c)/{d — c) < 20. Next, comparing the SIFs 
obtained with the Un and T„ representations in Table |^, we note that the results compare quite 
well, except when {d+c)/{d — c) ~ 1.0, and the discrepancy is bigger at the left-hand-side (LHS) 
than at the right-hand-side (RHS) crack tip. This occurs because of the "edge effect" [?]. If 42 
terms {i.e. A^ + 1 = 42) and T„ representation are considered for the case {d + c)/{d — c) = 1.01, 
then the normalized SIFs at the LHS and RHS crack tips are 3.6437 and 1.3302, respectively. 
Thus, when there is an "edge effect", the results are sensitive to the discretization adopted. 
Moreover, for the same number of collocation points, the level of accuracy attained with the 
Un representation is slightly different from that with the T„ representation. 

Figures ^ and ^ compare the crack profiles for f/„ and T„ representations. One may observe 
that the displacement profiles obtained from both representations practically agree within plot- 
ting accuracy, especially in Figure 0. Note that the displacement profile in Figure ^ is tilted to 
the left because of the "edge effect". Such effect is negligible in Figure |^. 
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Figure 5: The half plane of the antiplane shear problem for nonhomogeneous material with 
shear modulus G{x) = GqG^^. 

7.2 Mode III Crack Problem in Nonhomogeneous Materials [?] 

Consider the antiplane shear problem for the nonhomogeneous material shown in Figure |^ with 
shear modulus variation given by 



G{x) = Goef^- 



X 



(95) 



where Gq and f3 are material constants. Erdogan [?] has studied this problem in order to 
investigate the singular nature of the crack-tip stress field in bonded nonhomogeneous materials 
under antiplane shear loading. He uses a slope formulation^, while here we use a displacement 
formulation. To understand what can be gained through the displacement formulation, we first 
state the governing PDE and the boundary conditions for the crack problem: 



w(x,0) = , 



dx 



'yz 



[x, 0~^) = p{x) 



, — oo < X < oo 

X ^ [c, d] , 
X G (c, d) , 



y>o 



(96) 



^The governing integral equations are described by relationships (20), (21), and (22) of Reference [?] 
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where p{x) is the traction function along the crack surfaces {c,d); and because of symmetry, 
only the upper half plane y > is considered. By the Fourier transform we write w{x, y) as 

1 P°° 

w{x, y) = ^= [A{i)e^^^^y] e^^di , (97) 



where A{^) is to be determined by the boundary conditions, and 

im? = e + ^K . (98) 

Because of the far field boundary condition, \miy^oow{x^y) = 0, A(^) is found to have a non- 
negative real part which can be expressed as 



m = ^yVF+pe + e-^sm^Ve+W-e, m 



where the sign function S(-) is defined as 



1 , r]>0 
S{ri) = { , ri = (100) 

-1 , ri<0 . 



By applying the inverse Fourier transform to equation (pT]), one finds 



A(0 = -^ / w{t,0)e''^dt , (101) 

y/Zn Jc 



which leads to the following integral equation^]: 



G(x) f^ 
ay^{x,0)=p{x) = —-^ k{x,t)w{t,0)dt , (102) 



where 



/oo 
[A(Oe^(^)^] e*(*-")«de • (103) 

The trade-off between a displacement-based and a slope-based formulation can be seen here 
if one recalls the issue regarding decomposition of the K{E,) described in (|^). In summary, 

K{0 = lim ^e^(^)^ = -^ , if slope formulation is used, (104) 



and 



K{^) = lim A(^)e'*'*^^''^ = A(^) , if displacement formulation is adopted. (105) 



^Note that equations ( |102| ) and (103) correspond to equations (20) and (21) in reference [?], respectively. 
However, the present notation is different from [?], in which the dummy variable used for the Fourier transform 
is a, (a, 6) stands for the crack surfaces, and m{a) corresponds to our A(^). 

27 



The decomposition of K{$^) in (|104|) is difficult because of the term ^ in the denominator. On the 
other hand, the decomposition of K{^) in (|105|) can be achieved through a simple asymptotic 
analysis: 



Real part of A(0 = -^^V^^fW+i' '^^ -1^1 , 



(106) 



z X Imaginary part of A(0 = -^Sm\f7^^^W - e '^^ -y y ■ (107) 

Thus the governing hypersingular integral equation is found to be 



G{x) f^l 2 , /5 



27r J^ L(^ ~ ^y t — X 
where we have let 



+ -^— + A^(x,t) 



D{t)dt = p{x) , c<x <d , (108) 



D{t) = w{t,0) , (109) 



and the nonsingular kernel is 

N{x,t) 



r \ ~f3'VCcos[{t-x)C] 



-P^sm[{t-x)^] 



(p + v2SiP)^^/^^TWe-e) {^e +p'+ 2^/e+We 



d^ . (110) 



Recall that the function S(-) is defined by equation ( |100| ). As a consistency check, note that if 
/? = 0, then both the Cauchy singular kernel f3/{x — t) and the nonsingular kernel N{x,t) will 
be dropped from equation ( |10^ ) so that equation ( P5[ ) is recovered. 

Figure^ shows numerical results for displacement profiles considering a crack with uniformly 
applied shear tractions ayz{x, 0) = —po (|x| < a), and various values of the material parameter 
(3. Note that the cracks are tilted to the right because of material nonhomogeneity. Further 
numerical results, including SIFs at both tips of the crack and corresponding displacement 
profiles, are given by Chan et. al. [?] . From a numerical point of view, they have shown that 
essentially the same results are obtained either by f/„ or r„ representations for this specific 
problem [?]. 

7.3 Gradient Elasticity Applied to Mode III Cracks [?] 

One of the most relevant aspects of the formulas derived in Sections 4 and 5 is the evaluation 
of hypersingular integrals such as Ia{Tn,m,r) for m > 2 in the weight function W{s) given by 
equations (|10|) and ([T^). This example illustrates this point for the case m = 2. 
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Figure 6: The half plane of the antiplane shear problem for nonhomogeneous material with shear 
modulus G{x) = GqG^'^. The cracks are tilted to the right because of material nonhomogeneity. 



Paulino et. al. [?] have presented a hypersingular integral equation formulation for a mode 
III crack in a material described by constitutive equations of gradient elasticity with both 
volumetric and surface energy gradient dependent terms. A similar study, using a different 
approach, has been conducted by Vardoulakis et. al. [?] and Exadaktylos et. al. [?]. For this 
problem, the governing PDE is 



-r v*w + v'w = , 



(111) 



where i is the characteristic length of the material associated to volumetric strain-gradient 
terms, and w is the antiplane shear displacement. The boundary conditions are 



<Jyz{x,0) =p{x) 

IJ'yyz{x,0) = , 

w{x,0) = , 



\x\ < a 

— oo < a; < oo 
Ixl > a , 



(112) 



where the notation of references [?, ?, ?] is adopted here. Enforcing the governing equation 
( |111| ), imposing the boundary conditions ( [112| ), taking account of symmetry along the x-axis, 
and using Fourier transform method, Paulino et. al. [?] have obtained the following governing 
Fredholm hypersingular integral equation 



G 



n 



-2e 



1- 



'f_\2 
^21' 



{t - xY 



t — X 



Gi' 
+ k{x, t) )■ (j){t) dt H — —(j)\x) = p{x), \x\ < a (US) 
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with singularity a = 3, where the slope function 0(x) is defined to be 

, , , dw{x,0) ,^^ ^, 

<^(^) = Q^ (114) 

which satisfies the single-valuedness condition 

(p{t)dt = , (115) 



for the solution of the fracture mechanics problem. Here, k{x, t) is the nonsingular kernel given 
by 



k{x, t) = I ^ ^ . ^ 1 X ^ sin[e(t - x)]di ■ (116) 

where (—a, a) stands for the crack surfaces; I' is the characteristic length responsible for 
surface strain-gradient terms; G is the shear modulus of the material; p{x) is the known loading 
function; t is the integration variable, and x is the collocation variable. 

The behavior of the solution in terms of the density function 0(t), can be expressed as 

0(t) = i?(t)(a2-t2)| ^ (117) 

where R{t) is taken to be an expansion of Tchebyshev polynomials of first kind (r„). This 
example motivates the whole work of this paper because the analytical evaluation of the hyper- 
singular integral J3(T„,2,r) is needed for successfully solving the governing integral equation 
( |113|) . The expression for J3(T„,2,r) is given by (|55D . 

Note that the unknown density function is taken to be the first derivative of the displacement 
function, as described by equation ( p. 14| ) . For this particular example, the decomposition of the 
original kernel into singular and nonsingular parts, stated by equation (^), can be accomplished 
by means of partial fractions [?]. In general, this step of decomposition is not an easy task, 
as discussed in Example 2. An alternative is to consider the displacement function w(x, 0) as 
the unknown density function. In this case, a hypersingular integral equation with a = 4 is 
obtained. With order of singularity a = 4, the behavior of the density function D{t) in equation 
(I) can be expressed by [?] 

D{t) = R{t){a^-t^)l . (118) 

Thus one needs to evaluate the hypersingular integral /4(T„, 3,r) in order to implement the 
numerical approximation. The expression for /4(T„, 3, r) is given by equation 



The numerical results are presented in terms of SIFs (Table ^ and displacement profiles 
(Figure ^ by considering the slope function (f){x) as unknown in equation (|113|) and the T„ 
expansion to R{t) in equation (|117|) . Table |^ shows the convergence of the SIFs by choosing 
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Table 3: Stress intensity factors Kjjj{a) = 3 



N 



ifGEZ 



I a„, where a„ 



are the coefficients 



xpant 


on to R{t) in equat 


ion (ITT7) 


. Here £' 


= 






N 


£=0.8 


£=0.5 


£=0.2 


£=0.1 


£=0.05 


£=0.01 


£=0.005 


11 


20.3131 


15.8292 


7.4396 


4.5116 


2.6342 


0.1282 


0.0319 


21 


11.8757 


9.5632 


4.4791 


2.1538 


0.9541 


0.0898 


0.1602 


31 


11.6607 


9.3937 


4.3878 


2.0856 


0.9204 


0.1649 


0.0404 


41 


11.6665 


9.3983 


4.3902 


2.0878 


0.9246 


0.1378 


0.0682 


51 


11.6667 


9.3984 


4.3903 


2.0878 


0.9247 


0.1399 


0.0658 


61 


11.6667 


9.3984 


4.3903 


2.0878 


0.9247 


0.1400 


0.0653 


71 


11.6667 


9.3984 


4.3903 


2.0878 


0.9247 


0.1399 


0.0654 


81 


11.6667 


9.3984 


4.3903 


2.0878 


0.9247 


0.1399 


0.0654 


91 


11.6667 


9.3984 


4.3903 


2.0878 


0.9247 


0.1399 


0.0654 


101 


11.6667 


9.3984 


4.3903 


2.0878 


0.9247 


0.1399 


0.0654 



different values of volumetric gradient dependent term, £, and letting surface energy gradient 
dependent term £' = 0. The displacement profile in Figure |^ shows the "cusping" phenomenon 
at the crack tips, which is also described by Barenblatt's "cohesive zone" theory [?]. The 
interesting point is that the cusping crack obtained here is a natural outcome of the higher 
order gradient elasticity theory. Further numerical results and discussions are provided in 
reference [?]. 

8 Concluding Remarks 

Closed form analytical solutions are provided here for a broad class of improper integrals with 
hypersingular kernels and density functions approximated by means of Tchebyshev polynomi- 
als. Whenever possible, the symbolical and numerical tools of the computer algebra software 
MAPLE [?, ?, ?] have been used to verify the proposed solutions. A systematic approach for 
evaluating integrals when higher order singularities is also given in the present paper. 

The examples involve crack problems and aspects such as LEFM, nonhomogeneous mate- 
rials, and gradient elasticity theory (see Table [I|). All these problems are solved by means of 
Fredholm hypersingular integral equation formulations. When classical elasticity is used, both 
Tn and Un representations lead to essentially the same numerical results. For a crack problem 
in nonhomogeneous material, the difficulty that arises in splitting the singular and nonsingular 
parts from the original kernels can be circumvented by means of displacement-based, rather 
than slope-based, formulation. 

As a closing remark, we note that as material property variation in space and higher order 
gradient continuum theories are considered, the formulation of the crack problem and the 
associated kernels become quite involved. Thus, better analytical and numerical techniques 
are needed to successfully solve the governing singular integral equations. This paper is a 
combination in this sense. 



^Maple can evaluate some CPV integrals. However, in general, computer algebra systems are very limited 
with respect to hypersingular integrals. 
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Figure 7: Full crack displacement profile at d/a = 0.2 and i'/a = 0.1 under uniform crack 
surface antiplane loading ayz{x,0) = — (Tq; wq = acro/Go. 
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Appendix 

A Integrals associated with lower order n 

The general formula given in the text, e.g. equations (|2|), (|3|), (||), (§2]), ([S^), (§T]), (§g), and 
(ffOp are only valid above certain values of n. Thus the goal of this Appendix is to provide the 
expressions for integrals associated with lower order n. 



A.l /i(r,,l,r), n = 0, 1 

1 f' To{s)VT^^ 



IT j_i s — r 



ds = -r , \r\ < 1 (119) 



1 f^ Ti(s)VT^^ , . 1 , , , , 



A.2 Ii{T„,2,r), n = 0---5 



1 r W(l-.')»,^^^3_3^, |,|^i ,121) 



TT /_i s — r 2 



1 [^ Ti s l-s2 2 4 3 3,,, ,_^, 

_ / ^wv ^^s = /__^2^_ r <1 (122) 

TT Ji s-r 2 8''' ^ ' 



1 /-^ T2(s)(l 5')^.^^^2r5_4^3^9^^ 1^1^^ ^^23) 



TT /_i s — r 



TT Jli s-r 8 ' ' ' ^ ' 



/ ^ii^ ^^^ _ g^7 _ 20^5 + I6r3 - 4r , Irl < 1 (125) 



/ ^iil ^ds = IGr^ - 44r6 + 41r^ - 14r2 + 1 , Irl < 1 (126) 

J-i s-r 
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A.3 /i(r^,3,r), n = 0---5 

1 f To{s)il-s' 



2 15 5 o r, 

-as = r H — r — r 

TT I, s-r 8 2 



1 



(127) 



1 f T,{s){l-s^)l^ 5 15 , , 5 , , 



TT 



-1 



s — r 16 8 2 



:i28) 



1 f T,{s)il-s')l^_ 5^ 25^3 , . . 



TT 



-(is = — r r + 6r — 2r' 

s-r 12 4 



(129) 



1 /"i T3 s l-s2 2 25 55 2 ,,4 ,06 .8 

- -/- ^ ^^ —ds = \ r^ - 15r^ + 13r^ - 4r^ 

7r/_i s-r 32 8 



(130) 



1 / T4(s)(l s^)2 ^^ ^ _65^ ^ 20r3 _ SGr^ + 28r^ - 8r9 
n J_i s — r 16 



(131) 



1 f' T5(s)(l S^)2 ^^ ^ 31 _ ^^^2 ^ 55^4 _ 35^6 ^ go^8 _ ^g^lO 



TT 



s — r 
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(132) 



A.4 /i(t/„,3,r), n = 0---5 



1 /-^ f/os l-s2 2 15 535 

- f ^ ^^ —ds = r + -r^ - r^ 

TT / 1 s-r 8 2 



(133) 



1 f' [/i(s)(l-s2)l,_5 15 , , . . 



TT 



-rfs = r + 5r - 2r'^ 

s-r 8 4 



(134) 



1 f^ f/2(s)(l-s2)t ,25 3 5 , 7 

_ / ^wv ^^s= _^_ 10^3 + 11^5- 4r7 

TT J_i s — r 8 



(135) 



1 [^ f/3(s)(l-s2)f ,15 9 4 fi « 

_ 1 ^li^ ^(/s = + lOr^ - 25r^ + 24r^ - 8r® 

TT ./_i S-r 16 



(136) 



1 [^ f/4(s)(l-S 



TT 



2N^ 



-rfs = -5r + 30r-'^ - 61r^ + 52r^ - 16r^ 



s — r 



(137) 



1 f f/5(s)(l-s2)l 



n 



ds = l- 20r^ + 85r^ - 146r^ + 112r*^ - 32r^° 



s — r 



:i38) 
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A.5 /2(r^,3,r), n = 0---6 



1 [^ To(s)(l-s2)2 4 15 5 15 , , , , 

-'- °^ ^^ ' -ds = -5r^ + —r^-— , |r| < 1 (139) 



7ry_i (s-r)2 2 8 



1/ Ti(s)(l s2)2_^^^_g^5^^^^3_l^^^ 1^1^^ ^^^^^ 



TT y_;^ (s — r)2 4 



li^ T2(s)(l-s2)§^^^_^^^,^g^^,_75^2^5^ 1^1^^ ^^^^^ 



TT ./_! (s - r)2 4 2' 



1/ ZM£KL_^rfs = _32r^ + 78r5 - 60r3 + -r , |r| < 1 (142) 



1/ T4(S)(1 s')^^^^_72^8^igg^6_^80^4^gQ^2_65 ^ |^|^^ ^^^^^ 



TT y_i (s — r)2 16 

1 /-i T,{s){l-s 



TX 



f T is)i\ — sis 

=F ^V^ TT^ds = -160r^ + 480r^ - 510r^ + 220r^ - 30r , Irl < 1 (144) 

J^i [s-ry 



-1 ^i£li ^Ji^/s = _352rio + 1152r^ - 1386r^ + 730r^ - 150r^ + 6 , |r| < 1 (145) 

TT y_i {s-ry 
A.6 /2(t/n,3,r), n = 0---6 



1/ t/o(^)(l-^^)t,„_ ,^4 , 15^2 15 



, , ,^ -ds = -5r^ H r^ , Irl < 1 (146) 

TT ./_i (s — r)^ 2 8 



-/ ^^#^i^|^d5 = -12r^ + 20r^ - -r , |r| < 1 (147) 

vry_i [s-ry 2 



1 i' f/2(s)(l s2)2^^^_2g^g^^^^,_^^^2^25 



TT /_i Is — r)2 ' ' 8 ' 



- =L -llJl —^ds = -64r^ + 144r^ - lOOr^ + 20r , Irl < 1 (149) 

vr y_i [s-ry 
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- =L -^IJI ^^(/s = -144^8 + 364r^ - 305r^ + gOr^ - 5 , r <1 150 

7T J_i [s-ry 

-1 ^^^/^ ~ y' ds = -320r^ + 896r^ - 876r^ + 340r^ - 40r , Irl < 1 (151) 

vr J_i (s-r)^ 

-4 ^6(s)(l -s )^ ^g ^ -704r^o + 2160r^ - 2408r^ + 1155/ - 210r^ + 7 , |r| < 1 (152) 
n J_i {s-ry 

A.7 h{Tn,2,r), n = 0---4 

1 /-^ To(s)(l-s 



7r 7_]^ (s — r)' 



-ds = 3r , |r| < 1 (153) 



1 r' Ti(s)(l-s 



2N^ 



, , „ -ds = 6r^ , |r| < 1 (154) 

n .1^1 [s — r)-^ 2 



1 f^ T2(s)(l-s2)f 

' ^ ^^ ^ -rfs = 20r^ - 12r , |r| < 1 (155) 



TT y_i (s — r) 

1 r' T3(s)(i-s 



TT 






TT 7_i (s — r) 

A.8 /3(r„3,r), n = 0---7 



1 r^ T fsKl — s^^i 

4 ^V^ — ^rfs = 168r^ - 200r3 + 48r , Irl < 1 (157) 

J-i (s-ry 



If ^0(.)(l-.^)i^^^_^Q^3+15^^ 1^1^^ (^5g) 



TT y_;^ (s — ry 2 



1 /-^ Ti(s)(l-s2)f 4 2 15 , , . . 

-*- ^ ^^ ^ -rfs = -15r^ + 15r^ , |r| < 1 (159) 



TT J_i (s — ry 



1 [^ T2(s)(l-s2)| , o 75 , , , , 

-*■ ^ ^^ ^ rfs = -42r^ + 60r'^ - — r , |r| < 1 (160) 



TT 7_;^ (s — r)3 4 
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1/ ZM£KL^!)lrfs = _ii2r6 + 195r^-90r2 + f , |r| < 1 (161) 

1 i 4(g)(l-g )^ ^g ^ -288r^ + 588r^ - 360r^ + 60r , |r| < 1 (162) 
n J_i {s-r-y^ 

- 1 5(g)(l-s )^ ^^ ^ _720r^ + 1680r6 - 1275r^ + SSOr^ - 15 , Irl < 1 (163) 

-1 , ~ ^ ^ ' rfs = -1760r'^ + 4608r^ - 4158r^ + 1460r^ - 150r , |r| < 1 (164) 

TT y_i {s-ry 



TV J_i {s — r)^ 

-4224r^° + 12240r^ - 12768r^ + 5655r^ - 930r^ + 27 , |r| < 1 (165) 

A.9 /3(/7„,3,r), n = 0---6 

1/ f^o(.)(l-^^)^.rf, = _io.3 + l^,, |,|<1 (166) 



7r J_i {s — r)3 2 

1 i' ^^iMl^rf, = _30r^ + 30r^ - 1^ , |r| < 1 (167) 

71 J_i [s — r)-^ 4 

ds = -84r^ + llOr^ - 30r , |r| < 1 (168) 



1/ f/2(s)(l-.2)i,_ ,,^5 , ,,n.3 



TT y_i (s — ry 
1 r^ u (s)(i — s^)t 

- =L _l\^ —^ds = -22Ar^ + 360r^ - 150r2 + 10 , Irl < 1 (169) 

7rJ_i {s-ry 

1 f f/4(s)(l-s2)f 



TT 



T -^^ TT^ds = -576r^ + 1092r^ - 610r3 + 90r , |r| < 1 (170) 

J-i [s-ry 



1 _[ U5{s){l S )2^_^_ TAAn^S , QiQc^e oinn^4 , p;in^2 



TT y_i (s — r)' 



-rfs = -1440r^ + 3136r^ - 2190r^ + 510r^ - 20 , |r| < 1 (171) 



- =L ^iJl —^ds = -3520r^ + 8640r^ - 7224r^ + 2310r^ - 210r , Irl < 1 (172) 

n J_i (s - ry 
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A. 10 /4(r„,3,r), n = 0---8 



ds=-10r^ + -^ |r| < 1 (173) 



TT J-i (s-r)^ 2 ' 

1 [^ Ti(s)(l-s2)i 

- =^ ^1/^ TT^f^s = -20r^ + lOr , r < 1 174 

7rJ_i (s-r)4 



1/ '^^(')(^-'')lds = -70r' + 60r'-^, |r| < 1 (175) 



TT y_i (s - r)4 4 



1 r' T3(s)(l-s 



TT 7_i (s — ry 



ds = -224r^ + 260r'^ - 60r , |r| < 1 (176) 



-=L -^iJl ^^(is = _672r^ + 980r^ - 360r2 + 20 , Irl < 1 (177) 

7rJ_i [s-ry 



_[ 5[s){ ^ ''i(;s = _i920r^ + 3360r^-1700r3 + 220r , Irl < 1 (178) 



TT J_i {s — ry 

r^ T fsVi — s'^)^ 

T ^V^ — ^rfs = -5280r^ + 10752r^ - 6930r^ + 1460r2 - 50 , Irl < 1 (179) 

J-i [s-ry 



1 r^ T (s)(i — s"^)^ 

-=L ^±Jl — ^rfs = -14080r^ + 32640r^ - 25536r^ + 7540r^ - 620r , Irl < 1 (180) 

TT J_i {s-r-y 

I r^ T (s](l — s"^)^ 
-L ^^ '^ > rfs = -36608r^° + 95040r^ - 87360r'' + 33320r^ - 4560r2 + 104 , |r| < 1 

(181) 



TT J _i (s — ry 



A. 11 /4([/n,3,r), n = 0---6 

1 f' f/o(s)(l-s2)l 



rfs = -10r2 + -, |r| < 1 (182) 



TT J_i (s — r-y 

ds = -AOr^ + 20r , |r| < 1 (183) 



1/ U,{s){l-s^)l^_ ,^^3 



TT J_i {s — ry 
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1/ ^M£KL^!llrfs = _i40/ + ii0r2-l0, |r|<l (184) 



I f^ U (s)(l — s'^)^ 

' ^^ '^ ' -(18 = -448r^ + 480r^ - lOOr , |r| < 1 (185) 



TT J_^i (s — r-y 



1 r^ Ui{s){i-s 



f U^{s){l s )2 ^^ ^ _i344^6 ^ ^g2o/ _ eiOr^ + 30 , Irl < 1 (186) 

J -I {s-ry 



1 _[ f^5(s)(l S )^^„_ QO/in„7 I K070^5 onon^3 



vr 7_]^ (s — r)^ 



-rfs = -3840r^ + Q272r^ - 2920r'^ + 340r , |r| < 1 (187) 



r^ u is)i\ — s^^i 

=F ^V^ ^T^ds = -10560r® + 20160r'' - 12040r^ + 2310r2 - 70 , Irl < 1 (188) 

7-1 [s-ry 



A.12 Others 



1 ^1 o 



(l-r)5To(t)rft= / (l-r)5rft = -7r (189) 

1 J-i 8 



j\i - e)^T,{t)dt 



IT 



(190) 



(191) 



(192) 



Jjl - t')lut)dt = ^ 

/ (l-t2)5T„(t)rft = , for all n, and n^ 0,2,4. 
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